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ABSTRACT 

An analytical study is made of particle trajectories about a neutral sheet./ 
SuclT”a magnetically neutral sheet has been observed in the earth's magnetic? 
tail by Ness. 

A simple model is constructed where the magnetic field varies linearly, re- 
versing across the neutral sheet, and the electric field is everywhere homog- 
enous, perpendicular to the magnetic field and parallel to the neutral sheet. The 
important results in this case are that either protons or electrons never come 
out of the neutral sheet, and gain an unlimited amount of ener gy 

Another model is constructed by adding a small component of the magnetic 
field' perpendicular to the neutral sheet. Such a component could be furnished 
by the earth's dipole field in the tail. This component not only serves to bring 
the particles out of the neutral sheet, but turns both electrons and protons 90° 
from the accelerating electric field, in toward the earth . ! Furthermore, all par- 
ticles that are injected into the neutral sheet far enough from the earth emerge 
with small pitch angles to the magnetic lines of forceTj Particles injected at 
about 50 R e (earth-radii) emerge in a cone of about 1/2° about a line of force, 
the angle of the cone varying inversely as the cube of the distance from the 
earth. Protons injected at 50 R e emerge with an energy of about 13 kev, and 
electrons have only about 6 eV. In order to have emergent electrons with en- 
ergies of the order of auroral electron energies, the electrons must be injected 
at about 150 R e . Protons injected at this distance back in the tail will emerge 
with energies somewhere around 50 to 100 kev. 

When a magnetospheric magnetic field model is considered, consisting of a 
dipole field added to a tail field like the simple model, the above results imply 
that electron auroras will occur at higher latitudes than proton auroras and that 
proton auroras will occur before electron auroras in time. 

i Thes e results, showing that particles are accelerated in and then ejected 
from neutral sheets may also be important in situations such as solar flares, 
the day-side neutral sheet and neutral sheets occurring in interplanetary space, j 
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Cartesian co-ordinate system, e y is anti-solar, along the 
earth-sun line, e x points North, and e z is into the paper. 

The magnetic field reverses linearly across a neutral sheet 
of thickness 2d, the electric field is uniform and everywhere 
points out of the paper. R is the distance from the earth to 
the origin of the coordinate system 

Sketch of particle trajectories using the fields of the simple 
model (linear reversal of the magnetic field across the 
neutral sheet, no magnetic field component perpendicular 
to the sheet, and a uniform electric field in the -e z direction). 
Electrons are accelerated in the + e z direction, protons in the 
-e z direction. The amplitude of oscillation decays, so particles 
never come out and their energy goes to infinity 

Sketch of earth's dipole field lines. Near the geomagnetic 
equatorial plane (x = 0), the dipole field lines are predomi- 
nately in the + e x direction, and the magnitude falls off with 
the cube of the distance • 

Sketch of particle trajectories using the fields of the simple 
model plus a small component perpendicular to the neutral 
sheet. Such a component could be furnished by the earth's 
dipole field as indicated. Both protons and electrons oscillate 
about the sheet accelerating in opposite directions, and are 
turned in the same direction (toward the earth) by the small 
magnetic field component perpendicular to the sheet. When 
the particles are turned 90°, they are ejected from the neutral 
sheet. Electrons come out much sooner than protons, hence 
gain less energy. The dimensions shown are illustrative and 
not to scale 
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AURORAL PARTICLES FROM THE 
GEOMAGNETIC TAIL, I 


INTRODUCTION 

Ness (1964) has reported the discovery of a magnetically neutral sheet in 
the earth's geomagnetic tail, with measurements taken on board the IMP-1 
satellite, A somewhat uniform magnetic field of about 20 y (ly = 10 -5 gauss = 
10 -9 weber/m 2 ) is found from about 10 R e (earth-radii) to at least 30 R e in the 
anti-solar direction. This field is predominantly in the solar direction above a 
plane (roughly identified as the magnetic equatorial plane), and in the anti-solar 
direction below this plane. The magnetic field reverses across a sheet of thick- 
ness about 0.1 R e and goes to zero within this neutral sheet. 

It is of interest to look at charged particle trajectories about such a sheet, 
to see if the sheet can be a source of auroral particles, the particles originally 
coming from the solar wind. Adiabatic theory, which is good throughout the rest 
of the magnetosphere, cannot be used across such a neutral sheet, because the 
magnetic field changes significantly in distances much less than a gyroradius. 
The charged particle equations of motion must therefore be either solved ana- 
lytically, or computed numerically. 

Speiser (1964(b), 1965) has made numerical calculations, using a computer, 
of proton trajectories about a model neutral or current sheet. Using Liouville's 
theorem, he finds that the largest intensity of particles emergent from the sheet 
directly along magnetic lines of force, occurs in a thin output sheet. Such a 
model is thus capable of producing thin sheets of incoming particles for auroras. 
The mapping of these sheets onto the earth is near the auroral zones, and will be 
discussed in a future publication. 


A SIMPLE LINEAR MODEL 

The simplest model for the fields about a neutral sheet which can be dis- 
cussed analytically is: 


B = 



'y ’ 


( 1 ) 


1 


E = - a e 


(2) 


where b is the strength of the magnetic field when x = d, the sheet half-thickness, 
and a is the strength of the electric field. The electric field is added because it 
is believed that such a field exists across the magnetospheric tail. This belief is 
based on measurements of current systems and electric fields over the polar 
caps, and the assumption that such a field would be transmitted across the tail if 
the magnetic lines of force are equipotentials. 

The co-ordinate system being used and its relationship to the earth is 
sketched in Figure 1. The above magnetic and electric fields are also indicated 
in Figure 1. 



x 



Figure 1-Cartesian co-ordinate system, e y is anti-solar, along the earth-sun line, 
e x points North, and e z is into the paper. The magnetic field reverses linearly 
across a neutral sheet of thickness 2d, the electric field is uniform and every- 
where points out of the paper. R is the distance from the earth to the origin of the 
co-ordinate system. 


The equations of motion for a particle in the neutral sheet, using these 
fields are: 



x = Cj z x 

(3) 


o 

II 

: >> 

(4) 


z = ” C 3 - Cj x x 

(5) 

q b 

where C, = , and C„ 

1 m d 3 

q 

= — a. Mks units are used, 
m 
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Speiser (1964, (a),(b)) has given the solution to these equations for large 
time. (See Appendix A.) The result is that the particle executes a damped 
oscillation about x = 0 (the amplitude going as 1/ 1 1/4 ), while accelerating in the 
+ e z direction for electrons and the -e z direction for protons. This can be un- 
derstood as follows. The electric field in equation (5) accelerates a proton in 
the -e z direction, z then becomes proportional to -t ignoring the term -CjXx 
for the moment. Equation (3) then becomes 


X:-kx 


( 6 ) 


where k = Cj | z | , or k =^— j ^ t. 

Equation (6) is just the equation for oscillatory motion, with spring constant 
k. k gets larger with time, however, implying that the spring gets stiffer with 
time, thus the amplitude of oscillation decays. 

The term -CjXx in equation (5) is thus oscillatory and approaches zero, 
justifying its neglect, k is proportional to q 2 , so particles of either sign will 
oscillate. 

The oscillation in x(t) is due to the V x B force which is always toward x = o 
for x either positive or negative, because of the reversal of the magnetic field. 

The net result of this simple model is that charged particles of either sign 
never come out of the neutral sheet, and their energy increases without bound. 

See Figure 2 for a sketch of the results of this simple model. 


THE LINEAR MODEL WITH SMALL PERPENDICULAR 
FIELD ADDED 

Consider now the addition of a small magnetic field component perpendicu- 
lar to the neutral sheet, i.e., in the + e x direction referring to the co-ordinate 
system of Figure 1 . Such a component could be furnished by the addition of the 
earth's dipole field. (See Figure 3). The magnetic field of equation (1) now 
becomes: 


B = b 



x 

d 



(?) 


where t) will be assumed small. 
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Figure 2-Sketch of particle trajectories using the fields of the simple model (linear re- 
versal of the magnetic field across the neutral sheet, no magnetic field component perpen- 
dicular to the sheet, and a uniform electric field in the -e z direction). Electrons are ac- 
celerated in the +e z direction, protons in the -e^ direction. The amplitude of oscillation 
decays, so particles never come out and their energy goes to infinity. 
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Figure 3-Sketch of earth's dipole field lines. 
Near the geomagnetic equatorial plane (x - 0), 
the dipole field lines are predominantly in the 
+ e x direction, and the magnitude falls off with 
the cube of the distance. 


The equations of motion now become: 

x = Cj x z 

(8) 

y = c 2 77 z 

(9) 

z = -C 3 -CjXX~C 2 77y 

( 10 ) 


where Cj = (q/m) (b/d) , C 2 = (q/m) b, C 3 - (q/m)a,Cj and C 3 being the same as 
for the simple model. When 77 = 0 , the simple model results. 

Even without solving these equations, the particle motion can be understood 
qualitatively as was done for the simple model. 

Consider a proton (the arguments also hold for electrons with appropriate 
changes in sign) incident on this neutral sheet with small velocity. (Strictly 
speaking, the sheet is not now a neutral sheet.) The proton will be accelerated 
initially in the negative z direction (see Equation (10)), gaining a velocity, z , 
proportional to -t as for the simple model. From Equation (9), there will then 
be an acceleration in the -y direction proportional to z or -t, and thus y will 
be proportional to - 1 2 . (Note that y a C 2 C 3 a (q/m ) 2 , so either protons or elec- 
trons are turned toward the - e y direction.) 

As long as z is negative, Equation (8) will imply oscillatory motion in x(t) 
by the same arguments in the previous section on the simple model. Thus, as 
long as z is negative, the term -Cjxx in Equation (10) is oscillatory and will be 
assumed small as a first order approximation. The oscillations may now be 
either damped or growing depending on whether z increases or decreases with 
time. The last term in Equation (10) grows as + t 2 so that z will grow negatively 
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until z goes to zero, and z will then diminish in absolute value going to zero 
and even becoming positive after some time. Thus x(t) will execute damped 
oscillatory motion until z = 0 (the spring gets stiff er). After z becomes positive 
and until z = 0, x(t) will execute growing oscillatory motion. After z goes 
positive, however, x(t) will no longer oscillate but will increase exponentially, 
thus ejecting the particle from the neutral sheet. 

Summarizing, the addition of a small magnetic field component, + B x , to a 
linear neutral sheet model, serves to bring the particles out of the sheet and 
also turns charged particles of either sign in toward the earth (the - e y direction 
in Figure 1). Such a component +B X is consistent with the direction of the earth's 
field near the geomagnetic equatorial plane. Therefore, particles which have 
gained energy from the electric field are turned toward the earth, 90° from that 
field, and just as the energy gained is changed into motion toward the earth, the 
particle is ejected from the neutral sheet. It is shown in Appendix B that both 
protons and electrons emerge from the sheet with the same velocity. See Figure 
4 for a sketch of the particle trajectories in this model. 


RESULTS OF THE APPROXIMATE THEORY 

The detailed calculations are contained in Appendix B, to which the following 
alphabetic equation indices refer. First integrals of the equations of motion 
(Equations (8), (9) and (10)) are obtained exactly (Equations (i), (j) and (n)). A 
first approximation is used to obtain the time of ejection, r , which turns out to be 
inversely proportional to 77b(q/m) . (Equation (s)) Electrons are therefore 
ejected from the neutral sheet much sooner than protons. The velocity toward 
the earth at the time of ejection is found to be independent of (q/m), (Equation 
(w)) and the maximum pitch angle of the emergent particles (Equation (z)) is 
proportional to 77, and to |x 0 /u - c |, where u is the bulk flow velocity 
exterior to the neutral sheet, and c is a number which is not very different 
from 1. This result implies that a = 0 for 77 = 0, which means that the particles 
would never come out of the neutral sheet if the small perpendicular field were 
not added. This agrees with the result of the simple model, and implies that 
a will be small if 77 is small. Since a is also proportional to | x 0 /u — c | , par- 
ticles incident on the current sheet with x Q ~ u will all emerge with pitch angles 
close to zero. 


APPLICATIONS 

Choose R = 50 R e , b =20y, a = 0.3 volts/km. B x for the earth's dipole 
field atR = 50R e is B x % 0.3y. Therefore 77 = 1.5 x 10" 2 . Using Equation (v) for 
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Figure 4— Sketch of particle trajectories using the fields of the simple model plus a small 
component perpendicular to the neutral sheet. Such a component could be furnished by the 
4- .earth’s dipole field as indicated. Both protons and electrons oscillate about the sheet 

accelerating in opposite directions, and are turned in the same direction (toward the earth) 
by the small magnetic field component perpendicular to the sheet. When the particles are 
turned 90°, they are ejected from the neutral sheet. Electrons come out much sooner than 
1 protons, hence gain less energy. The dimensions shown are illustrative and not to scale. 
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the distance the particle has moved toward the earth at the time of ejection, we 
obtain y(r) ~ - 14 R e for protons and £ -40 km for electrons. The emergent 
velocity y (r) (Equation (w)) is about 1500 km/sec, giving an energy for protons 
of about 13 kev and about 6 eV for electrons. If the quantity | x /u - c | in 
Equation (z) is assumed to be of the order of 1, then in this example, the maxi- 
mum pitch angle is of the order of 0.3°. (All particles will emerge within a cone 
having a as the half-angle.) 

The right order of magnitude for auroral protons is 13 kev, but 6 eV is too 
small for auroral electrons. However, the output velocity is proportional to 
I/77, thus the output energy is proportional to I/77 2 . For the earth, 77 is pro- 
portional to 1/R 3 . So to get 6 kev electrons (which is of the order of auroral 
electron energies), we need only to have them injected in the neutral sheet at 
about 150 R e , assuming the neutral sheet extends that far. At that distance, the 
electrons still only stay in the sheet for about 6 R e before they are ejected. 

Injecting protons at 150 R e back in the tail would not increase their energy 
1000 fold as for the electron, because y(r) (Equation (v)) would need to be about 
-14,000 R e , so the protons would follow the neutral sheet in toward the earth 
until the stronger B x (larger 77 ) would bring them out. The proton energy would 
thus not get above 50 to 100 kev. 

If one adds the earth's dipole field to a tail field such as indicated in Fig- 
ure 1, then field lines from high latitudes are dragged back into the tail. A field 
line crossing the neutral sheet at 150 R e therefore comes from a higher latitude 
than a field line crossing at 40 R e . Therefore, the above results would imply 
electron auroras at higher latitudes than proton auroras, which is usually ob- 
served. Also, if some change in the characteristics of the solar wind (velocity, 
particle density, magnetic field direction) triggers this neutral sheet acceleration 
process, then the above results would predict proton auroras before electron 
auroras, since it will take time for the change in the solar wind to move from 
about 50 R e to 150 R e in the tail. The electrons emergent at 150 R e will all be 
within a cone of about .02° (Equation (z)). This implies that more electrons than 
protons should be able to get down to auroral altitudes. 
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APPENDIX A 


The first integrals of Equation (3) and (5) are: 


-^ c i ( x - x o> 


and 


(* 2 + ^ 2 ) + C 3 z = —• (i 2 + z 2 ) + C 


3 Z 0 


(a) 


(b) 


Equation (b) is just the equation of conservation of energy. The zero subscripted 
values refer to initial values. Equation (3) becomes, using (a): 

x = _ CjX^-z 0 +ic, (x 2 -x 2 ) +C 3 t) (c) 

For large enough time, the quantity in parenthesis on the right hand side of 
Equation (c) will be positive and monotonically increasing implying oscillatory, 
bounded motion in x(t). Thus for large time Equation (c) is approximately: 


x ~ - Cj C 3 x t 


-Gf 


a b 


x t 


(d) 


The solution to Equation (d) is given by Jahnke and Emde (1945, pp 147) as 


= /t 


' z v 3 (K /2 ) 


(e) 


where 



and Z v is a linear combination of Bessel Functions of the First and Second 
Kinds, of order one-third. Approximating (e) for large time, (Jahnke and Emde, 
1945, pp. 138), 


t-* 


A and B are constants depending on the initial values. 

z(t) may be obtained as a function of time by integrating Equation (a): 



Z(t) % - 



t + z r 


(g) 


+ smaller oscillatory terms 

From (f), after large time the amplitude of oscillation decays as l/t 1/4 , 
and from (g) and (b) it is seen that the kinetic energy increases as t 2 . 
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APPENDIX B 


First integrals of Equations (8), (9), and (10) are easily obtained. They are 



*2) 


= -C 3 (z - z 0 ) --(z 2 


~ z ' (j) - 


f y z dt 

Jo 


(h) 


y - y 0 = c 2 v ( z - z 0 ) 


(i) 


(z - i 0 ) = - c 3 t - -J (X 2 - Xq) - C 2 J] (y - y 0 ) 


( 3 ) 


Equation (j) becomes, on another integration: 


O - *o) = ‘ 


C 3 t 2 C 


- — J x 2 dt- c 2 1 ?J ydt + C 4 t 


X ai - | J ^ ^ t ^-4 

*'0 •'O 


(k) 


where 


C 4 = 


C 1 X 0 

+ -— ■ +C 2^y 0 


Using Equation (k) with Equation (i), we have: 
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i = y,t ’i 1* c s t ! + c 5 t - c, f x»dt-c|^[ ydt 

Jo Jq 


( 1 ) 


where 


C 2 C 3 C 1 C 2 

C - - - c - c c c = 

S ’ S - U 2 *“4 ’ '-'7 

2 2 


From Equation (9) we have: 


C 2 V f yz dt = i(y 2 -y 2 0 ) 
Jn 2 


(m) 


So that Equation (h) becomes the energy integral: 


I (i 2 + y 2 + i 2 ) + C 3 z =-L + y 2 + z 2 ) + C 3 z Q = 


(n) 


Using Equation (j), Equation (8) becomes: 


x = - k(t, x, y, rf) x 


(o) 


where 


C x 2 

k(t, x, y, tj) = C 3 t +—L — + C 2 77 y + C g 


(P) 


and 


C, x n 

C 8 = - — -C 2 V y 0 -i 0 . 

2 
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k goes to zero when 


t 




+ C 2 T) y + 



= r 


(q) 


t is therefore the critical time when k goes to zero and then becomes negative 
ejecting the particle. ( z becomes positive.) 


APPROXIMATIONS 

Integrating Equation (1), we obtain: 

• , C 6 t2 

y = y 0 + yo 1 +— — 


(r) 


Where the integral over x 2 in (1) is assumed small since x(t) is oscillating, 
and the integral over y in (1) is multiplied by rj 1 , so will also be neglected in 
this first approximation. 

Without loss of generality, but to facilitate finding the critical time of 
ejection r, initial conditions are chosen such that y 0 = y 0 = 0 and i Q = -CjX 2 /2, 
which implies that C. = C = C e = 0. Using Equations (q) and (r) with these 

4 6 8 

initial conditions, r becomes: 


From Equation (r) 



y(t) = - 


^c 5 t 3 

3 


(s) 


(t) 
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and 


y (t) = - v c s t 2 


(u) 


At t - t, Equations (t) and (u) become: 


y (T) = - 



and 


y (r) = - 


3a 

rjb 


(V) 


(w) 


Thus, in this first approximation, the ejection velocity is independent of 
(q/m) because electrons are ejected much sooner than protons. (As the next 
approximation, if y(t) from Equation (t) is used in the integral over y from 
Equation (k) it is found that r is increased by about 30%, y(r) is decreased by 
about 50%, and the ejection velocity is still independent of (q/m). The first 
approximation should therefore give the right order of magnitude.) 


PITCH ANGLE DISTRIBUTIONS OF EMERGENT PARTICLES 

In order to calculate the pitch angle of an emergent particle, we must esti- 
mate the maximum value x can have at the time of ejection. Since x(t) oscil- 
lates until t = r, it is reasonable to assume that x max ~ x 0 . Using this assump- 
tion and Equation (w) with Equation (7) for the magnetic field, with x = d, 
we have: 


V = k n e 

0 x 


3a ~ 

e 

rjb y 


(X) 


and the cosine of the pitch angle a , becomes: 





where u = a/b, the bulk flow velocity exterior to the neutral sheet, and only 
terms of order 77 2 are kept. For small a , cos a £ 1 - a 2 /2 , so 



The number 3 in Equation (z) comes from the 3 in Equation (w). This would 
be decreased to 3/2 by the second approximation mentioned above. 
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